Crossover from diffusive to strongly localized regime in two-dimensional systems 
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We have studied the conductance distribution function of two-dimensional disordered non- 
interacting systems in the crossover regime between the diffusive and the localized phases. The 
distribution is entirely determined by the mean conductance, {g), in agreement with the strong 
version of the single parameter scaling hypothesis. The distribution seems to change drastically at 
a critical value very close to one. For conductances larger than this critical value, the distribution 
is roughly gaussian, while for smaller values it resembles a log normal distribution. The two dis- 
tributions match at the critical point with an often appreciable change in behavior. This matching 
implies a jump in the first derivative of the distribution which does not seem to disappear as system 
size increases. We have also studied l/{g) corrections to the skewness to quantify the deviation of 
the distribution from a gaussian function in the diffusive regime. 
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I. INTRODUCTION 

The distribution function P{g) of the conductance of 
disordered systems is quite well understood in the metal- 
lic diffusive regime. It has a Gaussian shaped and so 
its first and second moments are sufficient to describe it. 
The second cumulant (irreducible moment) is a univer- 
sal constant-?' of order e^/ft- for a diffusive system of any 
dimensionality which explains the term universal conduc- 
tance fluctuations. 

On the other hand, the distribution function P{g) in 
the localized phase is known to be log normal in one- 
dimensional (ID) systems^, where the metallic regime 
does not exist, and recently has been found to be a 
Tracey-Widom distribution for the strongly localized 
phase in two-dimensional (2D) systems^. 

The crossover between the diffusive and localized 
regimes is possible both in quasi-lD and in 2D systems. 
Extensive studiea^>^ of quasi-lD systems have shown the 
distribution function to be independent of the system de- 
tails, with the average conductance being the only scal- 
ing parameter. This agrees with the one-parameter scal- 
ing hypothesis formulated for the mean conductance^ 
and later extended to the entire conductance distribu- 
tion function^. On the insulating side of the quasi-lD 
crossover, P{g) has been found to be essentially a "one- 
sided" log-normal distribution^ in agreement with ana- 
lytical studies of strictly ID systems.— 

The validity of the one-parameter scaling hypothesis 
for the conductance distribution function has also been 
studied in the localized phase of 2D systems^i^iiS where 
P{g) is now understood almost as well as in the diffu- 
sive phase. However, the crossover regime between the 
diffusive and the localized phases, where the localization 
length is of the order of the system size, is poorly under- 
stood in 2D systems. The numerical investigation of this 
regime is the focus of this paper. 

The conductance distribution can be described by its 
cumulants, which have been found in the diffusive regime 



to have the forro^ 



n < .go, (g) > 1 



(1) 



Here g is the dimensionless conductance (measured in 
units of e^/h), the brackets (...) stand for the ensem- 
ble averaging over disorder, (. . . )c denotes a cumulant 
and go is the mean conductance at the scale of the elas- 
tic mean free path i. It has been showni that there 
are no perturbative corrections to the second cumulant 
(variance). Therefore, in the absence of non-perturbative 
corrections it would remain universal (of order 1) for any 
value of {g). Since in the regime of strong localization, 
(g) <C 1, the variance should eventually decrease, and the 
non-perturbative corrections must exist. 

In the absence of a sharp transition between localized 
and delocalized states in 2D systems, it is possible that 
such non-perturbative contribution to the variance can 
be numerically traced already at the threshold between 
the diffusive and localized regimes where {g) ~ 1. As 
for the higher cumulants, in the metallic and crossover 
regime one can assess numerically only the third cumu- 
lant of the distribution, which represents the leading de- 
viation from the Gaussian. Equation Eq. ([1]) shows that 
it is proportional to 1/ {g) . The constant of proportional- 
ity has been calculated diagrammaticall}*ii and the result 
for quasi 2D systems is 
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-0.0020(.g)- 



(2) 



On the other hand, one can study numerically how 
the distribution function changes through the crossover 
regime from the almost Gaussian for (g) > 1 to al- 
most log- normal for {g) < 1. Analytical studies pre- 
dict also the appearance of the log-normal tails in the 
crossover regime^ which signify the emergence of pre- 
localizcd states^^. However, their statistical weight might 
be rather small in order to trace them numerically. 

Thus the aim of this paper is two-fold. First, we 
want to study systematically the full distribution func- 
tion P{g) in the crossover regime and check if the scaling 
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hypothesis extended to the full distribution applies in 
this region. Second, we want to obtain the corrections 
to the variance and the leading contribution to the third 
cumulant as a function oi l/{g) as we move away from 
the deep metallic regime into the crossover region. 

In the next section, we describe the model and the nu- 
merical procedure. In section III, we obtain the conduc- 
tance distribution function for a wide variaty of values 
of the disorder and the system size and present the re- 
sults on the applicability of single-parameter-scaling in 
the crossover region. In section IV, we calculate the sec- 
ond and third cumulants of the conductance distribution 
and analyze their \/{g) corrections. In the last section 
we summarize our findings. 

II. MODEL 

We have studied numerically the zero temperature con- 
ductance of the 2D Anderson model, described by the 
Hamiltonian 

i/ = ^ e^a\ai + t ^ ajoi -|- h.c. , (3) 

where the operator a\ (ai) creates (destroys) an elec- 
tron at site i of an square lattice and is the energy 
of this site chosen randomly between {—W/2, W/2) with 
uniform probability. The double sum runs over nearest 
neighbors. The hopping matrix element t is taken equal 
to —1, which set the energy scale, and the lattice constant 
equal to 1, setting the length scale. We have considered 
square samples of size Lx L. All calculations are done at 
an energy equal to —1, to avoid possible specific effects 
associated with the center of the band. 

The zero temperature conductance g is proportional to 
the transmission coefficient T between two semi-infinite 
leads attached at the opposite sides of the sample. 



where the factor of 2 comes from spin. From now on, we 
will measure the conductance in units of 2e^/h. We have 
calculated the transmission coefficient from the Green 
function, which was obtained propagating layer by layer 
with the recursive Green function methodic. This drasti- 
cally reduced the computational effort, and we can easily 
solve samples with lateral dimension up to 250 for the 
calculation of the distribution function, which requires a 
huge number of independent runs to get good statistics in 
the tails. The number of different realizations employed 
is of 10^ for most values of the parameters. We have 
considered wide leads with the same section as the sam- 
ples, which are represented by the same hamiltonian as 
the system, Eq. ([3]), but without diagonal disorder. The 
leads serve to obtain the conductivity from the transmis- 
sion formula in a way well controlled theoretically and 
close to the experimental situation. In the study of the 



distribution function, we use cyclic periodic boundary 
conditions in the direction perpendicular to the leads. 
In the calculation of the corrections to the variance and 
the skewness, we have also considered hard wall bound- 
ary conditions to make sure that they do not drastically 
change the results. The main conclusions are similar and 
we will present results for periodic boundary conditions, 
for which we have better statistics. 



III. CONDUCTANCE DISTRIBUTION 

We have obtained the conductance distribution func- 
tion for many values of the disorder and the system size, 
chosen in such a way that the system is in the crossover 
region, i.e., its mean value of the conductance is close to 
one. In figure 1 we represent many of these distributions 
as a function of the conductance for the values of the sys- 
tem size given in the legend. The disorder varies between 
W = 3 and 6. When the mean conductance is larger than 
one the distribution is basically a gaussian with an ap- 
proximately constant standard deviation, given by the 
value of the universal conductance fluctuations. As we 
approach the crossover, the shape of the high conduc- 
tance tail remains close to a gaussian function, but there 
is a drastic change of behavior of the distribution at the 
value g = 1. An important conclusion that can be ex- 
tracted from Fig. 1 is that the entire conductance distri- 
bution is uniquely determined by the mean conductance 
along the crossover region. This support the strong ver- 
sion of the single parameter scaling hypothesis. In several 
cases, we have adjusted the value of the disorder in such a 
way that the mean conductance is the same for different 
system sizes. In all these cases, the agreement between 
the different distribution functions is quite remarkable. 
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FIG. 1: (Color online) Conductance distribution function in 
the crossover between the diffusive and the localized regime 
regime for the values of L size indicated in the legend and for 
W ranging between 3 and 6. 
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For example, the distribution function for W = 4.2 and 
L = 60, whose mean conductance is equal 0.6677, is very 
similar to the distribution for W = 4.75 and L ~ 240, 
whose mean is 0.6686. 

For some distributions, we can appreciate a marked 
discontinuity in the first derivative of the distribution at 
g ~ I. This drastic change of behavior has been reported 
in different contexts, but it has never been analysed in 
depth for 2D systems. It appeared in the critical con- 
ductance distribution at the metal-insulator transition in 
3D systemai^, at the quantum Hall transitioni^, in the 
2D symplectic ensemble^^ and in quasi-one-dimensional 
systems^i^. In order to check if this apparent disconti- 
nuity in the fist derivative of P{g) is a finite size effect 
or not, we have run extensively three different systems 
with the same mean conductance, equal to 0.67 ± 0.01, 
and with sizes equal to 60, 120 and 240 (the values of the 
disorder are 4.75, 4.43 and 4.2 respectively). In Fig. 2 we 
plot the conductance distribution functions of these three 
systems. We first note the excellent agreement between 
the three of them, showing again that the distribution 
is enterily determined by a single parameter. Second, 
we appreciate that the possible discontinuity in the first 
derivative of the distribution does not seem to vanish 
away as the size increases. 
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FIG. 2: (Color online) Conductance distribution function for 
three different pairs of values of the disorder and the system 
size chosen so that the mean conductance is the same for all 
of them, equal to 0.67 ± 0.01. 



IV. 1/g CONTRIBUTIONS TO THE SKEWNESS 

The large amount of data gathered in our simulation 
allows us to study with enough precission how the con- 
ductance distribution deviates from a gaussian in the dif- 
fusive regime as a result oi 1/g corrections. We have 
quantify this deviation through the skewness, Sk, of the 
conductance. As we expect a leading contribution to Sk 




FIG. 3: (Color online) Skewness multiplied by (g) as a func- 
tion of l/{g) for periodic boundary conditions. Each symbol 
corresponds to a different system size. 



proportional to l/(g), we have plotted in Fig. 3 Sk{g) 
as a function of l/{g) for periodic boundary conditions. 
The data corresponding to the ballistic regime have not 
been drawn. The criteria to consider a data in the bal- 
listic regime was the ratio of system size to elastic mean 
free path, L/l, to be smaller than 20. We will discuss 
this criteria later on when analyzing the behavior of the 
variance. The elastic mean free path can be obtained 
through the dependence of the mean conductance with 
the system shape^i. For our model, we define I through 
the equation 



, , LJ 1 



(5) 



and Ly are the transversal and longitudinal 



where L. 

dimensions, respectively. We have fitted the data in Fig 
3 by an expression of the form 



Mi / \ b c 

{g) L 



(6) 



where the last term corresponds to a finite size contribu- 
tion. The straight line in Fig. 3 is the extrapolation to 
infinite size of Eq. ^ and intersects the vertical axis at 
—0.0005 ± 0.0004. This value is in reasonable agreement 
with the prediction of van Rossum et alr^ of —0.002 , 
if we take into account that this estimate does not con- 
sidered an specific set of boundary conditions and so we 
are not comparing exactly the same quantities. We note 
that in the region studied, the contribution to the skew- 
ness proportional to l/((?)^ is large and soon dominates 
over the negative linear contribution. 

We have also calculated the variance of the conduc- 
tance, = {g'^)c- In Fig. 4 we plot the variance of the 
data that we consider to be in the diffusive regime as a 
function of the inverse of the mean conductivity l/{g). 
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FIG. 4: (Color online) Variance as a function of l/{g) for 
periodic boundary conditions. Each symbol corresponds to a 
different system size. The horizontal line represents the value 
of the universal conductance fluctuations. 



The horizontal Une corresponds to the value of univer- 
sal conductance fluctuations for periodic boundary con- 
ditions, which is equaU^ (Tq = 0.1544. We note that the 
range of the vertical axis is small. There is a region in 
the ballistic regime where the variance presents a peak 
(not shown), whose magnitude depends on system size 
and is larger than the value of the universal conductance 
fluctuationsi^. In Figs. 3 and 4 we have only included 
data with L/l > 20, our criteria to be in the diffusive 
regime. We have calculated the eigenvalues of the trans- 
mission matrix, and for L/l < 20 there are many with 
transmission close to one. i.e., there are many ballistic 
channels contributing to the current. As g increases the 



variance gets closer to the universal value. 

We note in Fig. 4 that there is a systematic size de- 
pendence of the variance, breaking the strong version of 
single parameter scaling in this regime. This is proba- 
bly due to the ballistic regime, whose influence is much 
larger than expected. The existence of these flnite size 
effects together with the strong criteria needed to be in 
the diffusive regime, makes difficult a precise numerical 
determination of the universal conductance ffuctuations. 



V. DISCUSSION 

We have found numerically the behavior of the vari- 
ance and the skewness in the crossover regime is compat- 
ible with analytical predictions^. However, our results 
do not constitute a clear check since more precise numer- 
ical study of l/{g) corrections to the variance requires, 
on one hand, a very large number of realizations and on 
the other hand rather large sizes in order to avoid the 
ballistic regime. 

We have also observed the change in shape of the dis- 
tribution function when from diffusive through crossover 
regime. It can hardly be attributed, though, to the exis- 
tence of prelocalized statesii^^ which cannot be numeri- 
cally checked for available system sizes at present time. 
This implies that the pre-cxponential factor for these 
states should be relatively small. 
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